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Supersymmetric models provide many new complex phases which lead to CP violating 
effects in collider experiments. As an example, CP-sensitive triple product asymmetries 
in neutralino production e + e~ 



Xi Xi and subsequent leptonic two-body decays 



ini, £r — > XiA f° r ^ ~ e it J -i are studied within the Minimal Supersymmetric 



X? 

Standard Model. A full ILD detector simulation has been performed at a center of mass 
energy of yfs = 500 GeV, including the relevant Standard Model background processes, 
a realistic beam energy spectrum, beam backgrounds and a beam polarization of 80% 
and —60% for the electron and positron beams, respectively. Assuming an integrated 
luminosity of 500 fb^ 1 collected by the experiment and the performance of the current 
ILD detector, a relative measurement accuracy of 10% for the CP-sensitive asymmetry 
can be achieved in the chosen scenario. 



1 Introduction 

Supersymmetry (SUSY) [1] is among the most 
favoured and most studied extensions of the Stan- 
dard Model (SM) and is capable of solving many of 
its problems. One of its features is that the Mini- 
mal Supersymmetric Standard Model (MSSM) pro- 
vides a number of complex parameters which can 
serve as sources of CP violation. They are conven- 
tionally chosen to be the Higgsino mass parameter, 



/' 



\H\e 



the U(l) and SU(3) gaugino mass pa- 




Figure 1: Schematic picture of neu- 
tralino production and decay. 



rameters, Mi = \Mx\e^ and M 3 = \M 3 \e l ' p3 , respec- 
tively, and the trilincar scalar coupling parameters, 
Af = \Af\e^ A f , of the third generation sfermions 
(/ = b, t, t). CP phases can give rise to CP- violating 

signals in collider experiments [2] , which have to be measured to determine or constrain the 
phases independently of measurements of electric dipole moments (EDM). Although also 
CP-even observables, such as masses or branching ratios, are sensitive to the CP phases, 
CP-odd observables are needed for direct evidence of CP violation. 

In this report neutralino pair production e + e~ — > %° Xi, for i = 2,3, and the subsequent 
leptonic two-body decay of one of the neutralinos Xi ~^ ^ followed by Ir — > Xi ^ for 
t = e,n, at the ILC is studied [3]. Figure [1] shows a schematic picture of the process. The 
CP-sensitive spin correlations of the neutralino in its production process allow to probe the 
phase of the Higgsino mass parameter \i and the gaugino parameter M\ [4] . 
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A full ILD [5] detector simulation is performed in order to investigate in detail the 
prospects to measure CP-sensitive observables at the ILC. All relevant SM background 
is taken into account, simulated with a realistic beam energy spectrum and beam back- 
grounds [6j. 

2 CP-odd observables and benchmark scenario 

In neutralino production, effects from CP-violating phases can only occur if two different 
neutralinos are produced. CP asymmetries can then be defined with triple products of 
particle momenta. Due to the spin correlation the asymmetries show hints for CP phases 
already at tree level. For the process shown in Fig. [TJ a T-odd triple product of the beam 
and the final lepton momenta can be defined as [3] 

T = (p e - X P(+ ). (I) 

The corresponding asymmetry is 

_ (j(T > 0) - a(T < 0) 

AU) <j(T>0) + <r(T<oy [2) 

where a is the cross section for neutralino production and decay. Its sign depends on the 
charge of the leptons, which has to be tagged in the experimental analysis. 

For the full simulation study a benchmark scenario has been chosen such that the gaugino 
phase cf>i = 0.2-7T corresponds to a maximal CP asymmetry and the Higgsino phase is zero, 
since it is strongly constrained by EDM bounds. The other parameters in the neutralino 
sector are M 2 = 300 GeV, \M X \ = 150 GcV, = 165 GeV and tan/3 = 10. This leads to the 
neutralino masses m^o = 117 GeV. m^o = 169 GeV, = 181 GeV and m^o = 330 GcV. 

Xl 7 X.2 A3 A4 7 

while the slepton masses are mjr = 166 GeV and m| = 280 GeV. The neutralino pair 
production cross sections are calculated to be a(e + e~ — > X1X2) = 244 fb and er(e + e~ — > 
X1X3) = 243 fb, while the slepton pair production cross sections are a(e + e~~ — > e^e^) = 
304 fb and <r(e + e~ — > jj^jj^) = 97 fb. The slepton pair production is the main background, 
since there are two lightest neutralinos and two opposite-sign electrons or muons in the final 
state as in the case of the neutralino XiX? production. Furthermore, beam polarizations of 
(P e -,P e +) = (0.8,-0.6) have been chosen, which enhance slightly the SUSY cross section 
and the asymmetries, while the background from WW- and chargino-pair production is 
suppressed. In this scenario the CP asymmetries are A(T)^a^o = —9.2% and A{T)^>^o = 
7.7%. 



3 Detector simulation study and parameter fit 

The ILD is a concept under study for a multipurpose particle detector for the ILC. It is 
designed for an excellent precision in momentum and energy measurement over a large solid 
angle. A detailed description can be found in [5|. In the simulation all active elements and 
also cables, cooling systems, support structures and dead regions are taken into account [B]. 
The radiation hard beam calorimeter is used to suppress background from 77 events at low 
angles. All relevant SM backgrounds and SUSY processes are generated using Whizard [7]. 
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initial selection no significant activity in BCAL 
number of all tracks iVt rac ks < 7 

lepton selection £ + £~ pair with I = e, fj, 

| cos#| < 0.99, min. energy E > 3 GeV 
lower energetic £ with E < 18 GeV, or 
higher energetic £ with E > 38 GeV 
higher energetic ^ with S e [15, 150] GeV 

6> acop > 0.2tt, 6> acol > 0.2tt 

final preselection p™ iss > 20 GeV 
E vis < 150 GeV 
m« < 55 GeV 



Table 1: Preselection cuts, see Ref. [5] for details. 



3.1 Event selection and measured asymmetry 

A clean sample of signal events is needed in order to clearly measure the CP-violating 
effects in neutralino production. Otherwise the asymmetry will be reduced by the CP-even 
background events. Therefore, preselection cuts as listed in Tab. [1] are applied to reject as 
much background as possible, while preserving good signal efficiency. Electrons and muons 
are identified using the Particle Flow approach [3] . The cuts exploit the energy and angular 
distributions of the final state leptons, as well as the high missing transverse momentum 
p™ lss due to the escaping ncutralinos. Additional cuts on the total visible energy E v i s as well 
as on the invariant mass mu distributions further reduce the background contamination. 



Figure 2(a) shows the p™ lss distribution of the SM and SUSY background as well as 



of the signal after the lepton selection. It can be seen that most of the background is 



removed with the cut pJj? lss > 20 GeV. Figure 2(b) shows the distribution of the invariant 
di-lepton mass after all cuts except the one on m«. The signal lepton pair from X3 (X2) 
decays has a sharp endpoint at 51 GeV (22 GeV), which can also be exploited for mass 
measurements. The invariant mass cut also removes SM backgrounds from ZZ and WW 
production. The remaining event sample consists of 28039 X1X2 ^ XiXi^ 7^ r ) events, 
45966 X1X3 -> XiXi^ (* ^ r) events and 34223 U ->■ XiXi^ 7^ t) events. All other SM 
and SUSY background processes sum up to about 6000 events. 

In order to distinguish the X1X2 events from the X1X3 events, and to further clean the 
event sample, a kinematic selection procedure is applied, as described in [3j. A number of 
kinematic constraints derived from the final state momenta are used to classify events as 
signal or background. An event is selected only if it is classified exclusively as signal-like. 
Four event classes are considered: X1X2J X1X3J ^r^r an d W + W~ . Tablc[2]shows the number 
of events that are classified exclusively as one of the four event classes. It can be observed 
that the large contamination of the event sample by £r£~^ events can be drastically reduced. 

The CP asymmetry can now be calculated from Eq. <j2j) to be ^4(7^^o^o = — 11.3%±0.7% 
and A{T)y0^o = +10.9% ± 0.7%. The absolute values are slightly higher than the ones 
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Figure 2: (a) Missing transverse momentum p™ lss distribution of SM background, SUSY 
background and SUSY signal after the lepton selection, (b) Invariant mass mu distribution 
of the lepton pair after all preselection cuts except the cut on ma. The events are simulated 
for C = 500 fb -1 , beam polarization (P e - ,P e +) = (0.8, -0.6) at y/s = 500 GeV, and MSSM 
parameters as in the benchmark scenario discussed in Sec. O 



class 


only X1X2 


only X1X3 


only i + R t R 


only W+W~ 


X°iX°2 X?X?« (i ^ r) 


18343 


615 


51 


855 


x\x% -> X\X\U (£ + r) 


290 


20132 


372 


635 


all SUSY background 


1153 


3055 


5626 


951 


all SM background 


87 


256 


44 


81 



Table 2: Number of preselected events, that fulfill the requirements of the kinematic selection 
procedure, for C = 500 fb -1 . 



calculated in the benchmark scenario, since the asymmetry depends non-trivially on the cut 
values. This has been studied in [5] and can be taken into account in a parameter fit. 



3.2 Fit of the parameters in the neutralino sector 

In the final step of the analysis, the accuracy to determine the parameters in the neutralino 
sector of the MSSM is estimated. These are the six free parameters of the neutralino 
mass matrix \Mi\, M 2 , \fx\, tan/3, cf>i and As input for the fit a number of CP-even 
observables is used together with the measured asymmetries (see Ref . [3] for details) : m^o = 
117.3 ± 0.2 GeV, m^o = 168.5 ± 0.5 GeV, m^o = 180.8 ± 0.5 GeV, a(x°X°) x BR(xjf-> 
l R £) = 130.9 ± 1.4 fb, cr(x?X3) X BR(x° -> hi) = 155.7 ± 1.6 fb, a{xlxl) X BR(x° -> 
l R tf = 4.8 ± 0.3 fb, *{xlxi) x BR(xS] -> i R lf = 26.3 ± 0.7 fb and a(x° 2 X° 3 ) X BR(x° 2 -> 
i R €) x BR(x^ — > Zfl/) = 28.9 ± 0.7 fb. The fitted values of the parameters of the neutralino 
mass matrix are listed in Tab. [3] It is remarkable that the moduli of the phases <f>i, 0^ can 
also be determined with high precision, using the CP-even observables alone. However, only 
an inclusion of CP-odd asymmetries in the fit allows to resolve the sign ambiguities of the 
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I-Mil M 2 | M | tan/3 0! <^ 

150.0 ± 0.7 GeV 300±5GeV 165.0 ± 0.3 GeV 10.0 ± 1.6 0.63 ± 0.05 0.0 ±0.2 



Tabic 3: Results of the parameter fit. 

phases. Without the CP-odd asymmetries in the fit there is a twofold ambiguity, <f>\ = ±0.6, 
and even fourfold if 4>^ ^ 0. Thus, the triple product asymmetries are not only a direct test 
of CP violation, but are also essential to determine the correct values of the phases. 

4 Summary and conclusions 

The first full detector simulation study to measure SUSY CP phases at the ILC has been 
presented. Triple products of the final state lepton momenta in neutralino decays have 
been used as CP-odd observables. Realistic collider conditions have been simulated and 
all relevant SM backgrounds have been taken into account. A detailed cut flow analysis 
has been performed, including the development of a kinematic selection procedure that was 
used to obtain a very clean signal sample and to distinguish events from different neutralino 
decays. In the chosen benchmark scenario the asymmetry could be measured with a relative 
precision of 10% with 500 fb -1 of data. Finally, the parameters of the neutralino mixing 
matrix have been fitted to CP-even and CP-odd observables and the complex phases could 
be determined with a precision of about 10%. 
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